The g-factors of discrete levels in nanoparticles 
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Spin-orbit scattering suppresses Zeeman splitting of individual energy levels in small metal par- 
ticles. This suppression becomes significant when the spin-orbit scattering rate r^ 1 is comparable 
with the quantum level spacing S. The (/-factor exhibits mesoscopic fluctuations; at small St so it is 
distributed according to the Maxwell distribution. At St so — > the average (/-factor levels off at a 
small value g ~ (l/L) 1 ^ 2 given by the ratio of the electron mean free path I to the particle size L. 
On the contrary, in 2D quantum dots the (/-factor is strongly enhanced by spin-orbit coupling. 

PACS numbers: 71.24.+q, 73.23.Hk, 71.70.Ej 
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In bulk metals the spin-orbit interaction is known to 
affect the g-i actor of electrons. These corrections are rel- 
atively weak jjj ; even in such a heavy metal as gold, the 
(/-factor is 2.1, and differs only slightly from the nomi- 
nal value of 2. On the contrary, in recent experiments |^] 
with Al nanoparticles the measured values of the (/-factor 
were in the interval between 0.27 and 1.9. Later a similar 
effect was observed in gold nanoparticles ||, where the 
value g » 0.3 was measured. A systematic study of the 
effect of spin-orbit scattering on the g-factors of discrete 
levels was undertaken in Ref. [Q. In particular, it was 
demonstrated that when 4% of gold impurities are added 
to Al particle, the (/-factor drops from 2 to about 0.7. 

In bulk materials the (/-factor is conventionally deter- 
mined from the electron spin resonance (ESR) data. ESR 
involves transitions between states of the electron contin- 
uum. On the contrary, the experiments with nanoparti- 
cles |^-|| study the Zeeman splitting of individual elec- 
tron levels £,-: 
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Here H is the magnetic field, fiB = eh/2mc is Bohr mag- 
neton; m is the free electron mass. The splitting (Q) is 
linear in H as long as it is small compared to the quan- 
tum level spacing 6. The (/-factor determined by Eq. (Q), 
in general, varies from level to level. In the absence of 
spin-orbit interaction gi = 2. Evidently J§-|J] the spin- 
orbit scattering affects the values of gi. 

In this paper we develop a theory of the (/-factors of 
individual levels. Unlike its bulk value, the (/-factor de- 
fined by Eq. ([!]) is very sensitive to even weak spin-orbit 
interaction. The reason is that the magnitude of the cor- 
rection to g = 2 value measured by ESR in the bulk, is 
determined by the comparison of spin-orbit interaction to 
the Fermi energy Ep, whereas the relevant energy scale 
for (/-factors of individual levels is the quantum level spac- 
ing 5 <C Ep. Indeed, the spin-orbit interaction is usually 
described by the mean time of spin-orbit scattering r so . 
The scattering time should be compared with the time 



~ 1/6 that electron travels along a closed trajectory cor- 
responding to a quantum level. At <5r so 3> 1 the electron 
spin flips very infrequently, and the effect of the spin- 
orbit scattering is weak. On the contrary, at St so <C 1 
the spin flips the average of N — 1/St so times during 
the electron motion along the closed trajectory. Thus 
the average spin in such a quantum state is significantly 
less than 1/2, and the response (|l|) to the magnetic field 
is strongly suppressed. The (/-factor of this level can be 
estimated by assuming that electron spin-flips occur at 
random moments in time. The root-mean-square value 
of the spin is then 1/y/N, resulting in gi ~ ^/St so . 

The above estimate accounts only for the spin contri- 
bution to the (/-factor. In fact, the level splitting in a 
magnetic field is determined by the total magnetic mo- 
ment of the state 
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Here l z and s z are the operators of the angular momen- 
tum and spin, respectively; /i B = eh/2m*c, with m* be- 
ing the effective mass of the electrons. The brackets (■ • 
denote the expectation value in i-th eigenstate of an elec- 
tron confined to a nanoparticle in an infinitesimal mag- 
netic field in z-direction. In a nanoparticle of a generic 
shape, the orbital levels are not degenerate. Then in the 
absence of the spin-orbit interaction the time reversal 
symmetry dictates (l z )i — at H — > 0. It is impor- 
tant to note that in the presence of spin-orbit interaction 
('*)» 7^ 0, and the orbital motion contributes to the split- 
ting of levels by mag netic field f|. 

One can estimate the angular momentum of the elec- 
tron as the (directed) area A covered by its trajectory di- 
vided by the period of motion 1/8, that is, (l z )i ~ m*A8. 
To find A we notice that during time 1/Et the electron 
travels across the grain and, therefore, covers the area 
~ L 2 . (Here E? = D/L 2 is the Thouless energy, L is 
the grain size, and D is the diffusion constant for elec- 
trons in the grain.) During the period of motion 1/5 the 
electron bounces off the boundaries ~ Et/8 times. Since 
the direction of motion after each bounce is random, the 
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total directed area is A ~ L 2 ^JEt/5. Consequently, the 
root-mean-square angular momentum of each state is 



(3) 



Here I is the transport mean free path of electrons, and 
Uf is their Fermi wavevector. 

One can also see from Eq. (||) that in a three- 
dimensional (3D) nanoparticle the orbital contribution 
to the g-factor is small g° lh ~ y/l/L *C 2. However, in 
the case of strong spin-orbit scattering, <5t so <C 1, the 
above estimated spin contribution is also small, g^ p ~ 
yj St so <C 2. Therefore in this regime both contributions 
may have to be taken into account. 

We will now show that at i5r so <C 1 the g-factor of a 
level is a random quantity with the distribution function 
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where the averaging ((•••)) is performed either over an 
ensemble of 3D nanoparticles or over different levels in 
a single nanoparticle. Furthermore, we express ((g 2 )} in 
terms of two quantities, t so and /, which can be measured 
independently: 
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Here the dimensionless constant a is determined by the 
geometry of the nanoparticle; its exact value will be dis- 
cussed later. Eqs. (Q) and (||) are the main result of this 
paper. 

We begin the study of the g-factors of individual levels 
by finding relation between gi and the matrix elements of 
the operator of magnetic moment M. Due to the time- 
reversal symmetry, the levels in the nanoparticle are de- 
generate, with the wavefunctions and \Tipi), where 
T stands for the time reversal operator. In a small mag- 
netic field H the levels are split by the perturbation MH. 
Using the standard method of degenerate perturbation 
theory, one can easily find the splitting in the form (jl|) 
with the ^-factor 
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where the real vector jl is defined as 

M* + iliy = (Tiii\M\i>i), fi z = (il>i\M\ipi). (7) 
The distribution function p(0) is, by definition, 
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To perform the averaging in Eq. (ra), we use the Ran- 
dom Matrix Theory (RMT) approach H . Instead of the 



ensemble of nanoparticles with strong spin-orbit scatter- 
ing we will consider an ensemble of symplectic matrices 
of size 2N x 2N with N >■ 1. Then the eigenfunctions (ip\ 
and (Tip\ of the Hamiltonian are A^-component spinors: 
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The operator of magnetic moment M is a her- 
mitian matrix, which due to its time-reversal 
properties can be diagonalized to the form 
diag{Mi, -Mi, . . . , M k , -M k , . . . , M N , -M N }. In the 
basis of the eigenfunctions of M the matrix elements (Q) 
take the form: 
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(^|M|^)=^M fe (|^| 2 -| Xfe | 2 ), 

k=l 

N 

(T^|Af|V) =-2^M fc te- 
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The advantage of the RMT approach is that the ensem- 
ble averaging in Eq. (g) is easily performed using the 
Porter-Thomas distribution [HI of the matrix elements: 



A' 



n d 2 4> k d 2 Xk . 



2N(\4> k \ 2 + \ X kn 



(12) 



The averaging in Eq (g) with the help of (g_0|) — Q12|) re- 
duces to the calculation of N identical quadruple Gaus- 
sian integrals. The result has the form 
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In the limit of large N 3> 1, this integral becomes Gaus- 
sian too, and we find 
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Taking into account the relation (Jfy , we now immediately 
find the distribution function of the g-factor in the form 
(|J) with the mean square g defined as 
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As a phenomenological theory RMT enabled us to find 
the functional form of the distribution function; how- 
ever, the width of the distribution ( [To] ) is now expressed 
in terms of a phenomenological parameter TrM 2 /N 2 . 
The relation between this parameter and the microscopic 
properties of the system cannot be established within the 
RMT. To do this, one has to find an observable quantity 
Q which can be evaluated within both the phenomeno- 
logical RMT and a microscopic theory. 
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We choose Q to be the energy absorbed in unit time 
from a weak external ac magnetic field H(t) = Hq cosut. 
Such a field induces transitions between quantum levels 
in the grain, leading to the absorption of the energy in 
the grain. Within the RMT, the absorption can be found 
with the help of the Fermi golden rule as 
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The expression inside ((. . .)) is the rate of absorption of 
quanta of radiation interacting with the system. The 
absorption occurs due to the transitions from occupied 
states with k < kf to the empty states, p > kf, and are 
induced by the term ^H$Me~ lu)t in the corresponding- 
coupling Hamiltonian. 

In RMT the energy levels and the eigenfunctions are 
uncorrelated, i.e., the averaging over the energy levels 
and matrix elements in Eq. ( |16[) can be performed inde- 
pendently. Also, at uj ^S> S, one can neglect the correla- 
tions of the densities of states at energies separated by 
u), and we find 
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It is important to note that contrary to the above calcula- 
tion of the ^-factors of individual levels, Eqs. @ and (]?]), 
here the matrix elements M kp are between eigenstates 
with different energies. Upon averaging over the ensem- 
ble ({\Mk p \ 2 )) becomes independent of k and p. Its mag- 
nitude can be found by presenting the invariant TrM 2 as 
a sum of ((|Mfc p | 2 )) over all 2N values of k and p. Since 
the number of diagonal matrix elements 2 N is small com- 
pared to the number of the off-diagonal ones, AN 2 — 2N, 
we conclude that ((|M fcp | 2 )) = TrM 2 /4N 2 at iV » 1. 
Combining this relation with Eqs. (Jl?]) and (|l5|), we find 
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It is noteworthy that this result, obtained within the 
RMT, contains no phenomcnological parameters. It es- 
tablishes a relation between the property of a single level, 
the g-factor, and a macroscopic quantity Q insensitive to 
the effects of discreteness of levels. 

To evaluate the absorption rate Q for a given nanopar- 
ticle, it is convenient to express it as Q — ^ujA"Hq in 
terms of the imaginary part A" of the zz-component of 
the tensor of magnetic polarizability of the sample Aik, 
defined as Mi — AikHk, Ref. ||. Then the result for the 
mean square g-factor takes the form 
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It is well known || that at u> — > the imaginary part of 
the polarizability vanishes as A" {u>) cx uj. 

Since M — [i B l z + 2fi B s z , one can distinguish between 
the orbital and spin contributions to the magnetic polar- 
izability of the nanoparticle. The spin contribution has 
the form 



A a (w) = 
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Here the numerator is the usual static Pauli susceptibility 
of the electron gas of the nanoparticle, and the denom- 
inator accounts for the fact that spin correlations decay 
exponentially with the time constant r so 12. Substituting 
the imaginary part of the polarizability (go|) into (|l^) , we 
reproduce the first term of our main r esult ©. 

The orbital contribution to the magnetic polarizabil- 
ity is due to the magnetic moment of the eddy currents 
generated in the sample by the ac magnetic field. For a 
particle of a shape symmetrical with respect to the rota- 
tions around the z-axis one can easily find 
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Here a is the conductivity of the metal, V is the volume of 
the nanoparticle, and p\ is the "moment of inertia" of the 
grain, assuming unit density. For a spherical nanoparticle 
of radius L the combination of Eqs. ( |l9| ) and ( pil| ) repro- 
duces the second term in Eq. (|) with a = (6/5)(m/m*) 2 . 

For the thin ring geometry Eqs. ([l9]) and (^l|) result 
in ((g 2 )) = 3m 2 L 2 DS/n 3 h 3 , where L stands for the cir- 
cumference of the ring. In the context of the persistent 
current problem the orbital effect of magnetic field on 
the splitting of the energy levels was studied earlier by 
Kravtsov and Zirnbauer Q. They used the non-linear 
cr-model techniques to solve the general problem of 
crossover from the symplectic ensemble to the orthogo- 
nal one. A special limiting case of that solution gave rise 
to a distribution function of level splittings, which reas- 
suringly coincides with our result for the thin ring ge- 
ometry. Phenomenologically the crossover problem was 
solved within RMT by Mehta and Pandey S. Our ap- 
proach enables one to express their phenomenological 
crossover parameter to physical observables without re- 
sorting to cr-model calculations. 

In the case of weak spin-orbit interaction, 6t so ^> 1, 
the correction to the average ^-factor is small and can 
be found by perturbation theory. In the lowest order 
perturbation theory in spin-orbit coupling one finds W 
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2D quantum dot 



FIG. 



Ballistic nanoparticle 

Diffusive nanoparti cle 

1 ► 

1 M H/Stso 

1. Sketch of the dependence of the average g-factor 
on the strength of spin-orbit scattering, see text. We assumed 
m* jm ~ 1 in this plot. 

Assuming that the energy levels €j are equidistant, the 
sum ( p^) was evaluated by Kawabata flip] ]. However, in 
a disordered system the Wigner-Dyson statistics of the 
energy levels is a more realistic assumption. In this case 
one of the levels j can be close to level i, resulting in a 
particularly large correction to the g-factor. Taking into 
consideration the fact that the level repulsion in the or- 
thogonal ensemble suppresses the probability p of two 
levels being very close, p (ei — tj) = 7r 2 — e 7 - 1 /6<5 2 , M, 
we find that the average value of the sum in Eq. (P2|) is 
logarithmically large: 
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The logarithmic divergence in Eq 
energy scale 



I) was cut-off at the 
because of the additional 



level repulsion caused by the weak spin-orbit coupling. 

Our results are summarized in Fig. [IJ. The lower curve 
shows the dependence of the average (/-factor on the 
strength of the spin-orbit scattering in the case of a dif- 
fusive nanoparticle. In the regime of 1/i5t so <C 1 the 
behavior is described by Eq. (g3|) , and at 1/ 5t so 3> 1 the 
average g-factor drops in agreement with Eq. (|J) . In the 
latter regime ((g)) = (8/3-k) 1 / 2 ((g 2 )) 1 ^ 2 , as one can easily 
see from the distribution function (||). 

The values of the mean free path are usually not well 
known in the experiment, and it is possible that I ex- 
ceeds the grain size L. The middle curve in Fig. [j] shows 
schematically the behavior of a ballistic nanoparticle, 
where all the scattering of electrons is due to the reflec- 
tion from the boundaries. Qualitatively, in this case one 
can use the results for the diffusive case with the trans- 
port mean free path replaced by the size of the grain. It is 
then obvious from Eq. (||) that ((g)) levels off at a value of 
order unity in the limit of strong spin-orbit scattering, as 
illustrated in the figure. One can conduct a quantitative 
study in a simple model of a spherical ballistic nanopar- 
ticle of radius L with totally diffusive scattering off the 
boundaries. The magnetic polarizability can be calcu- 
lated with the help of the Kubo formula, which relates it 
to the correlation function of angular momenta at differ- 
ent times. This correlation function can be evaluated by 



considering the classical electron trajectories inside the 
nanoparticle, in the spirit of Ref. The exact result 

amounts to replacing I — > 5L/8 in Eq. (||) and keeping 
a the same as in the case of a diffusive sphere. Taking 
into consideration the rapid decrease of the distribution 
function (Q) of the g-factor at g 2 <C ((g 2 )), one should 
conclude that the nanoparticies showing the values of g 2 
well below the ballistic value ((g 2 )) = (3/4) (m/m*) 2 are 
most likely in the diffusive regime. 

The top curve in Fig. [TJ shows the behavior of the g- 
factor in the case of a 2D quantum dot in a perpendic- 
ular magnetic field. In accordance with Eq. (Q), in the 
strong spin-orbit scattering case, 1/St so ^ I, the orbital 
moment (l z ) reaches a very large value ~ \fk~F~l. In ex- 
periments with quantum dots in GaAs heterostructures 
the effect of the orbital field should be further enhanced 
due to a small effective mass to* ss 0.067m of electrons, 
so one can expect to find ((g)) ~ (m/m*)^/kpl. On the 
other hand, at weaker spin-orbit scattering, 1/5t so < 1, 
the orbital enhancement of the <?-factor is reduced, as the 
orbital motion does not affect the g-factor at all in the 
absence of spin-orbit coupling. The orbital contribution 
at 1/St so <C 1 can be found within the first-order per- 
turbation theory in the spin-orbit coupling, resulting in 
((g)) — (to/ 'm*)(hk F l j 'Stso) 1 I 2 . Thus, the g-factor of a 
2D quantum dot is dominated by the effects of the or- 
bital motion of electrons at 1/6t so > (kpl) -1 (to* /to) 2 , 
whereas at 1/5t so <C (kpl)^ 1 ^* /to) 2 values close to 
g = 2 should be observed. 
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